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ABSTRACT 

PSR J1048-5832 is a Vela-like (P=123.6 ms; r -20.3 kyr) 7-ray pulsar detected by Fermi, 
at a distance of ~2.7 kpc and with a rotational energy loss rate Esd ^ x 10^^ erg s~^. 
The PSR J1048— 5832 field has been observed with the Very Large Telescope (VLT) in the 
V and R bands. We used these data to determine the colour of the object detected closest 
to the Chandra position (Star D) and confirm that it is not associated with the pulsar. For 
the estimated extinction along the line of sight, inferred from a re-analysis of the Chandra 
and XMM-Newton spectra, the fluxes of Star D (F - 26.7; R - 25.8) imply a -0.13 < 
{V — R)o < 0.6. This means that the PSR Jl 04 8— 5832 spectrum would be unusually red 
compared to the Vela pulsar. Moreover, the ratio between the unabsorbed optical and X-ray 
flux of PSR Jl 04 8— 5832 would be much higher than for other young pulsars. Thus, we 
conclude that Star D is not the PSR J1048— 5832 counterpart. We compared the derived R 
and V-band upper limits (R >26A\ V >21 .6) with the extrapolation of the X and 7-ray 
spectra and constrained the pulsar spectrum at low-energies. In particular, the VLT upper 
limits suggest that the pulsar spectrum could be consistent with a single power-law, stretching 
from the 7-rays to the optical. 
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1 INTRODUCTION 

PSR J1048-5832 (B 1046-58) is a young, Vela-like pulsar in the 
Carina region, discovered during a 1.4 GHz Parkes survey of the 
Galactic plane (Johnston et al. 1992). It has a period P=123.7 ms 
and a period derivative P=9.63 x 10~^^ s s~^, which correspond 
to a characteristic age Tc ^^20. 3 kyr, a surface dipole magnetic field 
Bs ^3.5 X 10^^ G, and a spin down luminosity Esd ~2xl0^^ 
erg s~^. The NE2001 model (Cordes & Lazio 2002) yields a dis- 
tance of 2.7±0.35 kpc, based on the dispersion measure (DM). 
In the X-rays, PSR J1048-5832 was detected with the Rontgen 
Satellite (ROSAT) at 0.1-2.4 keV (Becker et al. 1997). Soon af- 
ter, observations with the Advanced Satellite for Cosmology and 
Astrophysics (ASCA) found possible evidence of extended X-ray 
emission (Pivoravoff et al. 2000), possibly associated with a pulsar 
wind nebula (PWN). High-resolution observations with the Chan- 
dra X-ray Observatory, indeed revealed the presence of an asym- 
metric PWN around the pulsar. No X-ray pulsed emission has been 



* Based on observations made with ESO Telescopes at the La Silla Paranal 
Observatory under programme ID 384.D-0386(A), 386.D-0585(A). 



detected so far by both Chandra and XMM-Newton, and a 3a up- 
per limit of 53% on the pulsed fraction between 0.5 and 10 keV 
has been determined (Gonzales et al. 2006). A recent re-analysis 
of the Chandra data of PSR J1048-5832 (MareUi 2012) could 
separate the pulsar emission from that of its PWN. See Marelli 
(2012) for details on the observations and data analysis. The com- 
bined Chandra plus XMM-Newton spectrum of the pulsar is fit- 
ted by a power-law (PL) with photon index rx=1.35di0.45 and 
A^H=(46.0 it 2.3) X 10^° cm~^, corresponding to an unabsorbed 
X-ray flux Px=0.490l°;3^2 x 10"^^ erg cm^^ s"^ in the 0.3- 
10 keV energy range. We note that, using the Chandra data, Gon- 
zales et al. (2006) tried to separate the pulsar and PWN compo- 
nents and derived a pulsar photon index Fx 2.4 (quoted without 
uncertainties) from the fit to the combined Chandra plus XMM- 
Newton spectrum. Our results differ from theirs since we used data 
from all the XMM-Newton detectors (both the EPIC-PN and MOS) 
and fitted the combined spectrum using the C- statistic approach in 
XSPEC. This allowed us to more efficiently separate the pulsar and 
PWN components in our spectral analysis and derive more accurate 
values for the pulsar photon index. Moreover, we used slightly dif- 
ferent radii to extract the pulsar counts: 20^^ instead of 25^^ for the 
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Table 1. Summary of the VLT R-band observations of PSR J1048-5832 
including: the exposure times (T), the number of exposures (N), the airmass 
sec(z) averaged over the N exposures, the seeing determined by the DIMM 
and its associated rms in (parentheses). 



Date 




N 


sec(z) 


seeing 


YYYY-MM-DD 


(s) 






n 


2010-12-04 


599.9 


5 


1.40 


0.80(0.02) 


2011-01-01 


599.9 


8 


1.47 


0.86(1.17) 


2011-01-04 


599.9 


5 


1.41 


0.91(0.09) 


2011-01-05 


599.9 


5 


1.25 


0.77(0.12) 



XMM-Newton data, and 2'' instead of V' for the Chandra data. PSR 
J1048-5832 was initially proposed by Kaspi et al. (2000) as the 
counterpart of the unidentified 7-ray source 3EG J 1048— 5840, de- 
tected above 400 MeV by the EGRET instrument aboard the Comp- 
ton Gamma Ray Observatory (CGRO). More recently, the Large 
Area Telescope (LAT) aboard the Fermi Gamma-ray Space Tele- 
scope detected for the first time 7-ray pulsations above 100 MeV 
from this pulsar (Abdo et al. 2009), highlighting a clear, double- 
peaked pulse profile very similar to the Vela pulsar. Deep optical 
observations of PSR J 1048— 5832 were performed with the Very 
Large Telescope (VLT) but only one object (Star D; V ^26.7) was 
detected close to the pulsar Chandra position (Mignani et al. 201 1). 
Here, we present new VLT observations of the PSR J 1048— 5832 
field in the R band that, together with the revised value of the ex- 
tinction along the line of sight obtained from our re-analysis of the 
Chandra and XMM-Newton data, confirm that Star D cannot be the 
pulsar optical counterpart. 



2 OBSERVATIONS AND DATA REDUCTION 

The R-band images of PSR J1048-5832 were obtained with the 
VLT Antu telescope at the ESO Paranal Observatory between De- 
cember 4, 2010 and January 5, 2011 (See Tab. 1) and are avail- 
able in the ESO archiv^B under programme ID 384.D-0386(A) and 
386.D-0585(A). All observations were performed in service mode 
with the FOcal Reducer/low dispersion Spectrograph (F0RS2, Ap- 
penzeller et al. 1998). F0RS2 is equipped with a red-sensitive MIT 
detector, a mosaic of two 2kx4k CCDs optimised for wavelengths 
longer than 6000 A. The pixel size of F0RS2 in standard reso- 
lution mode is 0^.^25 (2x2 binning), corresponding to a field of 
view of 8^3 x 8^3 over the CCD mosaic. Observations were per- 
formed in IMAGE mode, with standard low gain, normal readout 
(200 Kpix/s) and standard-resolution mode. For all the observa- 
tions, the target was positioned on CHIPl. The Rspec filter (A = 
6550 A; A A = 1650 A) was used for all observations. In order to 
minimise the saturation of bright stars and to allow for cosmic ray 
removal, a series of 23 short exposures of 599.9 s were obtained, for 
a total exposure of 13797.7 s (Tab. 1). Exposures were taken in dark 
time and photometric conditions. The average airmass was below 
1.5 and the seeing better than l^^ as determined by the Differential 
Image Motion Monitor (DIMM^ 

We used the standard packages in IRAF for the bias subtrac- 
tion and for flat field correction, selecting the closest-in-time bias 
and twilight flat field images available on the ESO archive. After 

1 www.eso.org/archive 

^ archive . es o . org/ asm/ambient- s erver 



these standard corrections, we aligned and stacked the reduced sci- 
ence images using the Swarp tool, included in the Scisoft 1.1 suit^E 
applying a 3cr filter to remove the hot/cold pixels and cosmic ray 
hits from the average image. Since all nights were photometric, we 
computed the photometry calibration of the average image by us- 
ing the average of the i?sPEC night zero points (< 0.05 magnitudes 
rms) and extinction coefficients computed over the four nights and 
the colour term for the September 2010-March 2011 semester. All 
these values are computed by the F0RS2 pipeline and are available 
through the F0RS2 data-quality control databas43- For the average 
science image we assumed the average airmass computed over all 
the exposures. We estimated that the effects of the airmass varia- 
tion over the different exposures (0. 1 1 rms) and the variation of the 
extinction coefficient (0.012 rms) over the four nights introduce an 
uncertainty of only 0.05 magnitudes on our photometry. 
The astrometric solution of the F0RS2 frame is based on the co- 
ordinates of the guide star used for the pointing. In order to im- 
prove the accuracy of the pulsar position on the frame, we re- 
computed the astrometric solution using the Two Micron All Sky 
Survey (2MASS) All-Sky Catalog of Point Sources (Skrutskie et 
al. 2006). Since the brightest stars in the field are saturated, we 
selected a subsample of 32 fainter, non-saturated reference stars, 
by excluding those at clearly bad position, e.g. close to the CCD 
edges. We then measured the star centroids through Gaussian fit- 
tingusing the Graphical Astronomy and Image Analysis (GAIA) 
tooO. We also used GAIA astrometric fitting routines to compute 
the pixel-to-sky coordinate transformation, which also account for 
the CCD distortions. The rms of the astrometric fits was ar 0Vl5 
in the radial direction. To this value we added in quadrature the un- 
certainty (Ttr=^' 1 of the registration of the F0RS2 image on the 
chosen astrometric reference frame. According to Lattanzi et al. 
(1997), this is crtr=yn/Nsa2MAss, where Ns=32 is the number 
of stars used to compute the astrometric solution, n=9 is the number 
of free parameters in the sky-to-image transformation model, and 
0-2MASS is the la mean positional error of 2MASS, which depends 
on the brightness of the reference stars. The selected 2MASS ref- 
erence stars have magnitudes 15.5 ^ K ^ 13, which corresponds 
to a2MASs ^0^.^2 (Skrutskie et al. 2006). We also accounted for 
the 0^.^015 accuracy of the link of the 2MASS coordinates to the 
International Celestial Reference Frame (Skrutskie et al. 2006). 
Thus, we estimated that the overall la positional uncertainty of our 
F0RS2 astrometry is ^0^.^18. For consistency, we also recomputed 
the astrometry calibration of the F0RS2 V-band image of Mignani 
et al. (2011) using the 2MASS catalog. The resulting rms of the 
astrometric fits was ar ^^0^.^16, leading to an overall la uncertainty 
on the V-band image of ^0''.^19. 



3 DATA ANALYSIS AND RESULTS 

As done in Mignani et al. (2011), also in this work, we assumed as 
a reference both the Chandra (a=10^48™12!64; ^=-58°32'03'.'60 
(MJD= 52859; Gonzalez et al. 2006) and radio interferometry 
(a=10^48™12^604; ^=-58°32'03'.'75 (MJD=50581; Stappers et 
al. 1999) positions that have radial positional uncertainties of 0^.^55 
and 0^.^08, respectively. These are shown in Fig[T] overlaid on the 
co-added F0RS2 R-band image (MJD=55566). Our astrometry is 



^ http://www.eso.org/sci/software/scisoft 
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Figure 1. VLT/FORS2 image of the PSR J1048-5832 field (R band; 
13797.7 s). The Chandra (Gonzalez et al. 2006) and radio interferometry 
(Stappers et al. 1999) positions are shown as circles. The star closest to the 
Chandra position is Star D of Mignani et al. (201 1). 

consistent with that of Mignani et al. (2011), although the accu- 
racy on our revised F0RS2 astrometry is better due to the better 
positional accuracy of 2MASS with respect to the GSC-2 (Lasker 
et al. 2008). As seen, there is no apparent candidate counterpart to 
PSR J1048— 5832 at the expected position. The object closest to 
the Chandra position is Star D of Mignani et al. (2011). We mea- 
sured its flux through standard aperture photometry, and obtained 
i?=25.8±0.2, where the associated error accounts for statistical er- 
rors and the accuracy of the absolute F0RS2 flux calibration (Mo- 
heler et al. 2010). No other star is detected at the Chandra position 
down to R ^26.4 (3a). We ignore whether Star D was the potential 
pulsar counterpart claimed by Solleman & Shibano\0 on the basis 
of the F0RS2 V-band data. Nonetheless, since Star D is also de- 
tected in the new F0RS2 R-band images, we used the whole data 
set to determine whether or not it can be a candidate counterpart to 
PSR J1048-5832. 



4 DISCUSSION 

Star D (V=26J±0.2) was ruled out by Mignani et al. (2011) as a 
possible counterpart to PSR J 1048— 5832 because of its ~ 3a off- 
set from the most recent Chandra position of Gonzalez et al. (2006) 
and an ^8% chance coincidence probability. However, since the 
time span between the Chandra (MJD=52859) and the F0RS2 V- 
band (MJD=55237) observations is ^6.5 years, the actual pulsar 
position could have been closer to that of Star D at the latter epoch, 
owing to its unknown proper motion. For instance, for a pulsar 
transverse velocity of 1200 km s~^, corresponding to a proper mo- 
tion of - 0'!9 yr"^ at a 2.7 kpc distance, the PSR J1048-5832 
position could have been closer to that of Star D by '-^0''.^6. Such 
a velocity is about three times as large as the average transverse 

^ http://archive.eso. org/wdb/wdb/eso/abstract/query?&progid=386.D- 
0585(A) 



velocity of radio pulsars (Hobbs et al. 2005). Pulsars with extreme 
velocities do exist, though, such as the Guitar Nebula pulsar, PSR 
B 2224-1-65, whose transverse velocity could be as high as 1600 km 
s"^ (Chatterjee & Cordes 2004). Although the probability that PSR 
J 1048— 5832 is moving at such an high velocity and exactly along 
a position angle of ^210°, i.e. that of the vector connecting the 
Chandra position to that of Star D, is small, we cannot rule it out 
a priori. If Star D was the pulsar counterpart, it should then have 
the same proper motion as hypothesised for the pulsar. We used 
the F0RS2 V and i?-band images, separated by 0.9 years, to 
measure a possible angular displacement of Star D through relative 
astrometry. We registered the two images using standard tools in 
IRAF, with an rms accuracy of ^ 0.05 pixel per coordinate. Un- 
fortunately, the faintness of Star D prevented us to fits its centroid 
with an accuracy better than 0.2 pixel per coordinate and we could 
only set a Icr upper limit of ^ 0.4 pixel on its displacement. Ac- 
counting for the F0RS2 pixel scale (0^25), this corresponds to an 
upper limit of ^ OVll yr~^ on the Star D proper motion. //"Star D 
was the pulsar counterpart, this would only allow us to rule out a 
transverse velocity > 1500 km s~^. 

We used the colour information on Star D, now available, to 
further investigate its possible association with PSR J1048— 5832. 
Star D has an observed {V — i?)=0.9=b0.3. We checked whether 
this colour would be compatible with a pulsar spectrum. We cor- 
rected the observed fluxes of Star D assuming the reddening to- 
wards PSR J 1048— 5832, computed from the hydrogen column 
density inferred from the X-ray spectral fits. The A^h derived from 
the best fit to the Chandra plus XMM-Newton X-ray spectrum 
(Marelli 2012) corresponds to a reddening E{B - y)=0.82±0.04 
according to the relation of Predehl & Schmitt (1995). This is much 
lower than the value inferred from the XMM-Newton spectrum 
alone (Marelli et al. 2011), E{B - V)=1.6toL and assumed by 
Mignani et al. (2011) in their analysis of the VLT data. Using the 
interstellar extinction coefficients of Fitzpatrick (1999), we derived 
Av=2.55±0A3 and Ak=1.91±0.09. Thus, if Star D were affected 
by the same extinction as the pulsar, its intrinsic colour would be 
(y— i?)o=0.26di0.33, i.e. its spectrum would be unconstrained. We 
note that the Vela pulsar has {V — R)o ^-03 (Mignani & Caraveo 
2001), where the uncertainty due to the reddening correction for 
A^H = (2.2zb0.5)xl0^° cm- 2 (Marelli etal. 2011) is marginal, i.e. 
bluer than it would be for PSR J1048-5832. However, Vela has a 
flat PL spectrum in the optical (Mignani et al. 2007), whereas some 
pulsars might have steeper spectra that rise towards longer wave- 
lengths, such as PSR B0540-69 (Mignani et al. 2010a), which 
would yield colours redder than Vela. Thus, the colour of Star D 
might not be incompatible with that of a pulsar. 

We note that the flux of Star D would be incompatible with the 
expected optical emission of PSR J1048— 5832 extrapolated from 
the luminosity of the Vela pulsar. PSR J1048-5832 has a spin- 
down luminosity Esd ^ 2 x 10^^ erg s"^, a factor of 3 smaller 
than the Vela pulsar (Esd ^ 6.9 x 10^^ erg s"^). Moreover, it 
is at a larger distance (2.7±0.35 kpc) than Vela (0.287=b0.02 kpc; 
Dodson et al. 2003) and is affected by a factor of ^-^20 larger red- 
dening. Assuming that the PSR J 1048— 5832 optical luminosity 
Lopt is the same fraction of the spin-down luminosity Esd as the 
Vela pulsar and accounting for the distance and reddening would 
then yield V ^32.2. If Star D were the PSR J1048-5832 counter- 
part, this would emit in the optical up to 3.7 x 10 ~^ of its spin- 
down luminosity. This value is ^100 larger than the Vela pulsar, 
and comparable to that of middle-aged and older pulsars (Zharikov 
et al. 2006). However, the evolution of the pulsar optical emission 
efficiency, r]opt = Lopt/ Esd, with the characteristic age is quite 
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Figure 2. Dereddened optical flux upper limits of PSR J1048— 5832 com- 
pared with the low-energy extrapolations of the X and 7-ray PLs that best 
fit the Chandra plus XMM-Newton (MareUi 2012) and FermifLAT (Abdo 
et al. 2009) data (thick black and red solid lines, respectively). The light- 
blue and yellow-shaded areas show the la uncertainty on the extrapola- 
tion of the X and 7-ray spectra, respectively. The uncertainty on the opti- 
cal flux upper limits due to the uncertainty on the reddening is ±0.04/iJy 
(±0.0004 X 10-27 erg cm-^ s-^ Hz-^). 



uncertain. In particular, this is true for spin-down ages between 10 
and 100 krs, where Vela is the only pulsar identified in the optical 
(Mignani 201 1). Thus, we do not know whether the sharp decrease 
of rjopt observed for Vela, which is only ^5 times older than the 
Crab but has a ?^1000 times lower efficiency, is peculiar of this ob- 
ject or is representative of all Vela- like pulsars. Indeed, the upper 
limits on the optical luminosity inferred for other Vela-like pulsars 
(e.g., Mignani et al. 1999; 2011; 2012) imply emission efficiencies 
riopt ^ 10~^, which would be consistent with that derived for PSR 
J 1048— 5832 //"Star D was its optical counterpart. 

We investigated this association from the corresponding X- 
ray-to-optical flux ratio Fx/Fy, where Fx and Fy are the un- 
absorbed Chandra and VLT fluxes of PSR J1048-5832 and Star 
D, respectively. For the assumed reddening, the association with 
Star D would imply, at most, an unabsorbed optical flux Fv ^ 
8.6 X 10~^^ erg cm"^ s"^ and Fx/Fy ^100. This value is 
much lower than expected from the comparison with other young 
rotation-powered pulsars (e.g., Zharikov et al. 2006), for which 
the non-thermal, unabsorbed X-ray-to-optical flux ratio is always 
larger than ^ 800. Thus, we conclude that Star D cannot be the 
PSR J 1048— 5832 counterpart. The R-band upper limit implies a 
flux F^ ^0-46 /iJy, corrected for the interstellar extinction. This 
flux value is quite close to that derived in the V band (F^ ^0-34 
/iJy). The optical flux upper limit on PSR J 1048— 5832 imply, for 
the revised value of the interstellar absorption, an optical emission 
efficiency r]opt ^ 1-9 x 10~^, accounting for the uncertainty on 
the A^H and the pulsar distance. 

We compared the optical flux upper limits with the low- 
energy extrapolation of the X-ray and 7-ray spectra (FigO. For 
the X-rays, we assumed the spectral model that fits the Chandra 
plus XMM-Newton data (Marelli 2012), a PL with photon index 
rx=1.35di0.45, while for the 7-rays we assumed the PL with pho- 
ton index r^=L38zb0.13 and exponential cut-off at ~2.3 GeV that 
fits the Fermi/LAT data (Abdo et al. 2009). As seen from FigH the 
extrapolation of the 7-ray PL can be consistent with the X-ray PL 



spectrum. This changes the conclusions reported in Mignani et al. 
(2011), based on the XMM-Newton spectrum, where the presence 
of a spectral break between the X-rays and the 7-rays was appar- 
ent. FiglSalso shows that the optical flux upper limits are just above 
the extrapolation of the X-ray and 7-ray spectra. This means that 
the multi- wavelength spectrum of PSR J1048— 5832 could be con- 
sistent with a single PL stretching from the 7-rays to the optical 
(see also Durant et al. 2011), at variance with other Fermi pulsars 
(e.g., Mignani et al. 2010b; 2011; 2012), where there is evidence 
for spectral breaks between the PL components in the different en- 
ergy bands, which are possibly representative of different energy 
and particle distributions. Thus, the detection of PSR J1048-5832 
in the optical would be crucial to confirm this scenario. 



5 CONCLUSIONS 

The new VLT observations of PSR J1048-5832 confirm that the 
object detected closest to the Chandra position (Star D) is an unre- 
lated field star. No other possible counterpart is detected down to 
R -26.4 and V -27.6. Thus, PSR J1048-5832 remains unidenti- 
fied in the optical. It is clear that its detection is close to the sensitiv- 
ity limit of 8m-class telescopes. Observations at longer wavelength 
might be more successful, though, were the optical/infrared spec- 
trum of PSR J1048-5832 follow the extrapolation of the X-ray 
PL. An updated Chandra or radio position for PSR J1048-5832 
would allow one to minimise the position uncertainty due to the 
unknown pulsar proper motion and would represent an advantage 
in the search for its optical/infrared counterpart. 
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